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In this paper we study the symmetry protected Majorana edge states for the Z2 topological order 
of the Wen-plaquette model and the toric-code model and calculate the dispersion of the Majorana 
edge states. For the system with translational symmetry, the Majorana edge states are gapless and 
have the nodal points at k = and k = ir. For the edge states of the toric-code model without 
translational symmetry, the edge modes become gapped. 
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I. INTRODUCTION 

Topological properties of topologically ordered states 
can be partially characterized by their gapless edge 
states pHj]. For example, the fractional quantum 
hall(FQH) states possess robust gapless edge modes as 
protected by the energy gap of the bulk, which can be 
derived from the Laughlin wave function [5|. On the other 
hand, people can obtain the bulk properties from the in- 
formation of the edge states by bulk-boundary correspon- 
dence. For example^dge excitation of FQH states is chi- 
ral Luttinger liquid [l|, [2( while that of non-Abelian FQH 
states @ is more exotic and related to (l+l)-dimensional 
conformal field theories d, Q- 

Recently, several exactly solvable spin models are 
found with non-Abelian topological ordered state or Z2 
topological ordered state as the ground states, such as 
the toric-code modelQ], the Wen-plaquette mo del pi. Hoj 
on a square lattice and the Kitaev model on a honeycomb 
latticeQ. In non-Abelian topological order of the Ki- 
taev model on a honeycomb lattice the elementary exci- 
tation becomes non-Abelian anyon with nontrivial statis- 
tics. And the edge state is ID gapless chiral Majorana 
modes. However, for a simpler example of a topological 
order - Z2 topological order, gapless edge states have not 
been found while a gapped edge state may exist [TTI - [l3 | 
instead. In the toric-code model, it is pointed out that 
there are two distinct types of boundaries (the smooth 
one and the rough one) with gapped edge states (Til 
We call them all zigzag boundary in this paper. 

In this paper, instead of considering the zigzag bound- 
aries that have gapped edge states, we study the Wen- 
plaquette model and the toric-code model with a new 
type of smooth boundary (which differs from the the 
smooth edge in Ref.flll HH) and use Majorana formu- 
lation to derive the effective theory of the gapless edge 
states. See the illustrations below. The gapless edge 
states are protected by translational symmetry along the 
boundary When the translational symmetry is broken, 
the edge modes will have energy gap. 

The paper is organized as follows. In Sec. II, we study 



the edge states of the Wen-plaquette model. In this sec- 
tion, we give the definition of the edge states by string op- 
erators and use the Majorana representation to describe 
Majorana edge states. In Sec. Ill we study the edge states 
of the toric-code model by similar approach. Finally, the 
conclusions are given in Sec.V. 



II. THE EDGE STATES OF THE 
WEN-PLAQUETTE MODEL 

A. The Wen-plaquette model 

The Wen-plaquette model is defined on square lattice 
with the Hamiltonian 
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and g > 0. af, a? are Pauli matrices on site i. 

The ground states of the Wen-plaquette model are 
known as Z 2 topological state [l,[loj]. The ground state is 
denoted by Fj = +1 at each plaquette. For this model, 
the elementary excitations are Z2 vortex (to type exci- 
tation denoted by Fj = — 1 at even sub-plaquette) and 
Z2 charge (e type excitation denoted by F, = — 1 at odd 
sub-plaqucttc). In addition, there is a mutual-scmion 
statistics between Z2 vortex and Z2 charge. A Z2 vor- 
tex and a Z2 charge annihilate with each other into a 
fcrmionic Z2 link-excitation which is a pair of Z2 vortex 
and Z2 charge. The fermions have flat band - the energy 
spectrum is F(k) = 4g, which implies that they cannot 
move at all. 

The ground states of Z2 topological state have topo- 
logical degeneracy. Under the periodic boundary condi- 
tion (on a torus) , the ground states of the Wen-plaquette 
model have four-fold degeneracy on even-by-even (e * e) 
lattice, two-fold degeneracy on even-by-odd, odd-by-even 
and odd- by-odd lattices. For a system on a cylinder, the 
ground states have two-fold degeneracy. Physically, the 
topological degeneracy arises from the presence or the 
absence of ir flux of fermions through the hole, as illus- 
trated in Fig. [TJ 
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FIG. 1: The illustration of a system on a cylinder. The topo- 
logical degeneracy arises from the presence or the absence of 
7r flux of fermions through the hole. 



B. String representation of the edge states for the 
Wen-plaquette model 

In Ref . [HI , it is pointed out that for the Wen-plaqucttc 
model with smooth open boundary condition, there ex- 
ist gapless edge states. Fig. 2 shows the Wen-plaquette 
model with a smooth open boundary and Fig. 3 shows 
the Wen-plaquette model with a zigzag open boundary. 
For a finite L x x L y lattice with a periodic boundary con- 
dition only along y-direction, there are two edges along 
y-dircction. This model can be realized by setting g = 
for one column of plaqucttcs and it remains exactly solv- 
able. The ground states have ~ 2 Ly -fo\d degeneracy and 
can be viewed as gapless edge excitations on both bound- 
aries described by Majorana fermion. These gapless edge 
states can be mapped to a Majorana fermion system with 
flat band exactly. Above argument comes from the ex- 
actly solvable model. But how the edge states change 
when there exist external fields? From the effective Z2E 
type mutual U(l)xU(l) Chern-Simons(CS) theories of 
Z2 topological order of the Wen-plaquette model, we have 
shown that there are right-moving and left-moving gap- 
less edge excitations described by Majorana fermions, 
provided that the edge is in the x- or ^-direction(l5j]. 
The presence of the translational symmetry in the x- or 
y-dircction is crucial for the existence of the gapless edge 
excitations for the Z2E type mutual U(l)xU(l) CS the- 
ory and the lattice model. In addition, from the classi- 
fication of the Z 2 topological order on a square lattice, 
we found that the emergent Majorana edge states of the 
Wen-plaquette model will always gapless at k x = 0, 7r (or 
k y =0, 7r) in momentum space 

To describe the edge states of the ground states (the 
planar codes), we define the fermion string operators. 
In the bulk, the fermion string operators are W/(C) = 
H &i m EU ) where C is a string connecting the middle of 



FIG. 2: The illustration of a smooth boundary of the Wen- 
plaquette model. The boundary has translational symmetry. 
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FIG. 3: The illustration of a zigzag boundary of the Wen- 
plaquette model. The boundary has no translational symme- 
try. 



the nearby sites, and i m are sites on the string. l m = z 
if the string does not turn at site i m . l m = x or y if 
the string makes a turn at site i m . l m = y if the turn 
forms a upper-right or lower-left corner. l m — x if the 
turn forms a lower-right or upper-left corner. Taking the 
open boundary condition into account, the loop C can 
be different. For system with open boundary condition, 
C is a string from one boundary to another. It is obvious 
that there are two kinds of strings, one has two ends at 
the same boundary (we call it C s ), the other terminates 
at different boundaries (we call it Cf). See Fig. 4. 
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FIG. 4: The illustration of three fermion strings: the closed 
string C 3 corresponds to the string with two ends at the same 
boundary; another closed string operator Ci corresponds to 
the string with two ends at the different boundaries; the third 
corresponds to the open fermion string with one end (the blue 
spot) on the boundary and the other (the green spot) in the 
bulk (site Io in the bulk). There exists Majorana fermion 
mode on each end of the open string. 



Now we can create the edge states |edge) by perform- 
ing a fermion string operation Wf(C s /i) connecting the 
boundaries on the ground state |0) as 

ledge) = W f (C s/l ) |0) . 

Both the edge states |edge) and the ground states |0) can 
be denoted by Fi = +1 at each plaquette and have the 
same ground state energy as 

E Q = H w |edge) = H w |0) . 

It is obvious that for the Wen-plaquette model, both 
types of fermion strings connecting the boundaries are 
condensed. These edge states have exact zero energy or 
flat band. 

Let's explain why the edge states are really Majorana 
fcrmions. People know that for an open fermion string, 
there exists a Majorana fermion modes at each end. In 
Fig. 4, we have shown an open fermion string, of which 
one end (the blue spot) is on the boundary, the other (the 
green spot) is in the bulk (site 1$ in the bulk). The end 
of fermion string in the bulk corresponds to a Majorana 
fermion mode. Since the fermion parity of system is con- 
served, there must exist another Majorana fermion mode 
at the boundary of the system which is another end of the 
fermion string (site i on the boundary). Thus each end of 
the fermion string connecting the boundaries corresponds 
to a Majorana fermion mode. On the other hand, each 
fermion string connecting the boundaries corresponds to 
a two-level state denoted by W f (C s/l ) |0) = ± |0) . Then 
we can say that each end of the fermion string connect- 
ing the boundaries is really a Majorana fermion mode. 



There are two types of fermion string operators, W/(C S ) 
denotes the string with two ends at the same boundary, 
Wf(Ci) denotes the string with two ends at the different 
boundaries. 

When we add external field terms 
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the degeneracy of the edge states will be removed and we 
will get the dispersive edge states. In the following parts 
we will calculate the dispersion of these edge states and 
point out that for the Wen-plaqucttc model with smooth 
boundary the gapless edge states are protected by the 
translational symmetry. 

In addition, for this case, the edge states on a fi- 
nite L x x L y lattice with a periodic boundary condi- 
tion only along y-direction always have two-fold de- 
generacy. The two-fold degeneracy is characterized by 
the closed string operator along .T-direction Wf(Ci) = 
af + , x <rf +2 - x • • • a* +L ^. For the case of W f (C{) |0) = |0) , 
there is no 7r-flux inside the hole and the Majorana 
fcrmions on the edge have periodic boundary condition; 
for the case of Wf(Ci) |0) = — |0) , there is a 7r-fiux in- 
side the hole and the Majorana fermions on the edge have 
anti-periodic boundary condition. 



Majorana edge states for the Wen-plaquette 
model 



In the following parts we will study the (symmetry 
protected) edge states by Majorana representation. Now 
we may use a Majorana mode 7* to denote an end of the 
fermion string at boundary. For the boundary shown in 
Fig. 2, the corresponding effective Hamiltonian of a single 
edge mode is given by 



edge 



where 7, is the Majorana operators at edge position i, 
and obeys algebra relation {7^,7^} = (7t)' = 7i- 



1 . Quantum tunneling effect of Majorana modes 



Firstly, we select the Majorana edge states character- 
ized by fermion string operator W/(C S ) = Yl <V"- A 

m 

fermion string operator Wf(C s ) that connects the two 
points on boundary can be considered as quantum tun- 
neling processes of virtual quasi-particles moving along 
the path. The quantum tunneling process of fermions is 
defined as : at first a single (bulk) fermion and an edge 
fermion are created together. Then this bulk fermion 
propagates and disappears at the boundary site j. And 
a string of ct™ is left on the tunneling path behind the 
virtual fermion, that is just a string operator Wf(C s ). 
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and get 




FIG. 5: The illustration of the relation between the hopping 
parameter Ji in the effective model of the edge states and the 
quantum tunneling process. There is translational symmetry 
for the edge states along the boundary. 
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where is the excited state of an edge fermion mode 
and a bulk fermion at link I (See Fig. 5). From H j^j) = 
(E Q + Ag) , we have 

( l —^H!)\l) = (-^-)|*<). 

E o -H ~V 

Then the bulk fermion move from /-link to / + 1-link by 
an operation, then turn around to / + 2-link by 

an , « operation. Finally, the bulk fermion disap- 

pears by performing the operation at I + 2-link by co- 
operation at site i + 1. 

Then we can get the energy shift of the state (1| as 



SE^J^iMHd—^Hjyil) 

j=o Eo - H 



h x (h*) h* 
(-4<7) 3 ' 



(7) 



For simplicity, we use |0) and |1) to describe the two 
degenerate eigenstates of the string operator Wf(C s ), as 



Wf(C a )\0) 



, W>(C.)|1) = -|1) 



With the perturbation term Hi, the quantum tunnel- 
ing processes occur - the fermion propagates from one 
end of the string to the other. Thus we can use 
the quantum tunneling theory in Ref.[l9|, [2(| to ob- 
tain the hopping parameters for the Majorana modes, 
Jij. The nearest neighbor hopping parameter J\ of 
the Majorana edge modes corresponds to the short- 
est fermion string with two ends at the same bound- 
ary, of which the fermion string operator is Wf(C s ) = 
a*af + ^&y +ix + - y a* + - y . See the illustration in Fig.5. 
And the next nearest neighbor hopping parameter J2 
corresponds to a fermion string with two ends at the 
same boundary, of which the fermion string operator is 

Let's calculate J\ by the higher order perturbation ap- 
proach. From the higher order perturbation approach in 
Rcf. 19, 20] , we obtain the energy shifts of two eigenstates 
|0) and |1). We take the energy shifts of the quantum 
state |1) as an example which is 



Using the same approach we can get the energy shift of 
the state (0| as 

(8) 

Finally an energy difference e of the two quantum states 
is obtained as 



e = SE 1 - SE 



2h x (h z ) 2 h y 
(-4 5 ) 3 



(9) 



which is the strength of the nearest neighbor hopping 
parameter J\ of the Majorana edge modes, 



2h x (h z ) h y 
(-4<?) 3 " 



Similarly we can derive the next nearest neighbor hop- 
ping parameter Ji of the Majorana edge modes by using 
the same approach as 



Ji = 



2h x (h z yhP 
(-4 5 ) 4 ' 



(10) 



2. Symmetry protected Majorana edge states 



6E 1 = (l\Hi( 
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-Hif^n) 



(5) 



To calculate ( 

starting point of the generation of an edge fermion mode 



En— Hi 



Hi) |1) , we can choose site i as the 



In this part, we use the Majorana representation to 
derive the dispersion of the Majorana edge states. By 
the following representation, 

7, = ^ (a + 4) , 
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3. Interference between the Majorana fermions on two 
boundaries 
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FIG. 6: The excited energies of a single edge mode of the 
Wen-plaquette model. The nodal points of Z2 topological 
order are fixed at k = and k = n (the red spots). We have 
set J 2 = -0.2 Ji. 



we have 

-Hedge — ^ ^ ^ 

= ~Y^2( Ci + 4)( c i+i + c l+i) 

+ % ~Y Y,^ + C i)(Ci+2 + 4+2) + - 
i 

%J\ ^—^ \ 1 + 4-4- 

i 

+ J2^ Cl + 2 ~ C i+2 Ci + 44+2 + C1C1+2] 

i 

+ ... 

In the momentum space we derive the energy spectrum 
of the edge states as 

E+ ~ 2(Jisinfc+ J 2 sin2fc), (12) 
= 0. 

Here the quantum states with zero energy are un- 
physical. Thus the excited energy of the edge states is 
given by 

AE = \E+\ = 2|Jisinfc + J 2 sin2fc| . 

We found that the nodal points of Z 2 topological order 
are fixed at k — and k = ir on an edge which is pro- 
tected by the Z 2 topological invariants and translational 
invariance. See Fig. 6. 



In this part we study the interference between the Ma- 
jorana fermions on two boundaries. Now the effective 
Hamiltonian of two coupling Majorana edge modes is 
given by 



-Hedge — i JtlA,ilB.i 
i 

+ 1 ^2 J\lA,ilA.i+l + i X] J llBSlB,i+l 

i i 

+ i ^2 J 2lA,ilAS+l + i ^2 J 2lBSlB,i+l ( 13 ) 

i i 



where A and B denote the indices of the two boundaries. 
Jt represents the coupling strength between two bound- 
aries. 

To derive J t , we need to consider another type of string 
operators as W f {Ci) = ^ +iy af +2iy ■ ■ ■ a z i+Lyiy . Such 
string operator is described by the quantum tunneling 
process of a link-fcrmion moving from one boundary to 
another along y-direction. The energy difference of the 
double states is 



e = -8g(—) L * -> J t . (14) 



Now L y is not very big. 

In addition, for a finite system we need to consider the 
quantum tunneling effect along cc-direction which will re- 
move the two-fold degeneracy of all quantum states in- 
cluding the edge modes. Because such quantum tunnel- 
ing effect is characterized by the closed string operator 
W f {Ci) = Sf+^CTf+a^ • • •Oi+L.e, alon S x-dircction, the 
energy splitting is about — &g{-z^) L * ■ Here we consider 
the case with big L x . Thus we can ignore this quantum 
tunneling effect. 

For this case of Wf(Ci) |0) = |0) , there is no 7r-fiux 
inside the hole and the Majorana edges have periodic 

boundary condition. From 7^ = ( c& A + v A A /V% and 
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1b, i = (c B ,i + c B ^j /y/2, we have 

-Hedge — * X! JtHASlBS 
i 

+ i ^2 JllA.il A.%+1 + i ^ J llB,ilB,i+l 
i i 

+ i 22, J 2lA,ilA,i+l + i 22 J 2lB,ilB,i+l 
i i 



i J i 
~2~ 

iJi 
~2~ 

iJ 2 
~2~ 

iJo 



J2 { CA >* + C A, 
i 

n 



n 



CB.i + C' B 



CAA + C 



J2 (''" ■ + c k 



CB,i + C B 
CAA+l + 



t ^ 

B,i+1 J 



CBA+1 + C 



CA,i+2 + C' Ai+2 



CBA+2 + C\ 



,«+2 



(15) 



In the momentum space we have 



-Hedge = ^2(c A ,k, c\_ k , C B ,k, C B _ k )e(k) 



( 4* \ 



k>0 



CA-k 
c B,k 

V c B ,-k J 



where 



i{k) 



e k e k iJt/2 iJ t /2 

e k e k iJt/2 iJt/2 

-iJt/2 -iJt/2 £k £k 

-iJt/2 -iJt/2 e k e k 



(16) 



(17) 



and = Ji sin k+ J 2 sin 2k. Now we have the eigenvalues 
as 



E = 0,0,2 ( Ji sin k + J 2 sin 2k) ± J t 



(18) 



In particular, for periodic boundary condition. The mo- 
mentum is 

2irn 

k = , n = 1.2..., M 

where M is an integer number. 

On the other hand, for the case of W(Q) |0) = - |0) , 
there is a 7r-flux inside the hole and the Majorana edges 
have anti-periodic boundary condition. Wc have similar 
energy levels. However due to the anti-periodic boundary 
condition, the momentum is 

2irn — ir „ 
k= M , n=\,2...,M 

where M is an integer number. The two cases are always 
degenerate for an infinite system. 
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FIG. 7: The excited energies of two edge modes of the Wen- 
plaquette model with interference effect. The nodal points of 
Z2 topological order shift away from k = and k = -k (the 
red spots). We have set J2 = — 0.2Ji, J t = 0.2J\. 




FIG. 8: The excited energies of two edge modes of the Wen- 
plaquette model with interference effect. There exists a gap- 
less point with quadratic dispersion (the red spot). We have 
set J 2 = -0.2 Ji, Jt = 1.08 Ji. 



Because the quantum states with zero energy are un 
physical, the excited energies are 

E + = |2(Jisinfc 
£L = |2(JisinA: 



- J2 sin2fc) 

- J2 sin2fc) 



■Jt\, 
-Jt\. 



(19) 



We found that the nodal points of Z 2 order are not fixed 
at k = 0/ir on an edge due to the interference effect be- 
tween the Majorana fermions on two boundaries. Sec 
Fig. 7. For J t > I.O8J1, the edge states become full 
gapped by mixing the edge states on different bound- 
aries. See Fig. 9. 

In particular, for a special case Hi = h z ^2&f, the 

z 

bulk fcrmion can only move straightforwardly but can- 
not turn a corner. Now we have J\ = J 2 = 0. Then 
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FIG. 9: The excited energies of two edge modes of the Wen- 
plaquette model with interference effect. The edge states are 
full gapped. We have set J2 = — 0.2Ji, Jt = 2Ji. 



the effective Hamiltonian of two coupling Majorana edge 
modes becomes 



H e d g e = iJt / ] lA,ilB,j- 



(20) 
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FIG. 10: The illustration of a smooth boundary of the toric- 
code model. The boundary has translational symmetry for 
the case of A = B. When there is no translational symmetry 
for the case of A 7^ B, the unit-cell along the edge has two 
sites. Here we denote the unit-cells by i and i + 



By the following representation of Majorana modes, 

lA,i = (<* + 4) /V2, 72M = ( c * - 4) /(V2»). 
we have 

-Hedge = iJt X!( Ci + C i)( _ *)( C i _ 4) 

i 

This is a complex fermion system with two energy levels. 



III. THE EDGE STATES OF THE TORIC-CODE 
MODEL 

The toric-codc model is described by the Hamiltonian 



Htc = -A J2 % - B J2 %i+ ( 21 ) 

zGeven iGodd 



where 



>i CT i+e x "i+e^ +e y " i+e v ■ 



fr x >y> are p au Ji matrices on sites, i. In this paper we only 
consider the case of A > 0, B > 0. The ground state of 
the toric-code model is also Z2 topological state which 



is similar to that of the Wen-plaquette model. For the 
toric-code model, if we choose A > 0, B > 0, the ground 
state of it is denoted by Z, t = +1 and X { = +1 at each 
plaquette. 



String representation of the edge states for the 
toric-code model 



To describe the edge states of the ground states (the 
planar codes), we also define the fermion string opera- 
tors. In the bulk, the fermion string operators are also 
Wf(C) = Y[ ^i™ i where C is a string, and i m are sites on 

m 

the string. l m = y if the string does not turn at site i m . 
l m = x or z if the string makes a turn at site i m . l m = x, 
if the string turns around the X-plaquette; l m = z, if the 
string turns around the Z-plaquette. For system with 
open boundary condition, C is a string from one bound- 
ary to another. It is obvious that the two types of string 
operators correspond to the string with two ends at the 
same boundary (we call it C s ) and that with two ends 
at different boundaries (we call it Ci), respectively. Sec 
Fig. 5. Similarly we can create the edge states |edge) by 
doing a fermion string operation Wf {C s n) connecting the 
boundaries on the ground state |0) as 



ledge) = W f (C s/l ) |0) . 
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B. Majorana edge states for the toric-code model 




Now we may use a Majorana mode 7,; to denote an 
end of the fermion string at boundary for the toric-code 
model. For the boundary shown in Fig. 10, the boundary 
has translational symmetry for the case of A = B, When 
there is no translational symmetry for the case of A 7^ 
B, the unit-cell along the edge has two sites. Here we 
denote the unit-cells by i and i + 1. Now we have two- 
sublattice in a unit-cell and denote them by 1 and 2. 
The corresponding effective Hamiltonian of a single edge 
mode is given by 



-Hedge = Ul,a X] 7l,»72,i + iJ\,b X! 72,»7l,t+l 
i i 

+ 1J2 ^2 7i,t7M+i + iJ 2 72,i72,i+i- (22) 



FIG. 11: The illustration of the relation between hopping 
parameters Ji, a , Ji.b in the effective model of the edge states 
for the toric-code model and quantum tunneling processes. 
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FIG. 12: The excited energies of a single edge mode of the 
toric-code model. The edge states are full gapped without 
translational symmetry. We have set J2 = — 0.0125Ji, n and 
Ji,b = 0.5Ji, a . The spectrum splitting comes from a small 
next nearest neighbor hopping. 



The nearest neighbor hopping parameters J\ a and Jij, 
of the Majorana edge modes correspond to the fermion 
string operators W f (C s ) = ofof^ ofj^ j_ s of, - and 



W f (C s ) = a?al +s J? + ^ + ,a v i+i 



respectively. See the 
illustration in Fig.ll. And the next nearest neighbor 
hopping parameter J2 corresponds the fermion string op- 
erator W f (C s ) = ^+eJl + e x +e^i+^ + 2^1 + 2e y ' Simi " 
larly we can derive the nearest neighbor hopping parame- 
ter Ji and the next nearest neighbor hopping parameter 
J2 of the Majorana edge modes by using the same ap- 
proach in above section as 



2h x {h y fh x 2h z {h y fh z 



(-2A-2B) 3 ' 1,0 (-2A-2B) 3 
2h x (h y fh z 
2 ~ (-2A-2B) 4 ' 

See Fig.ll. 

By defining 71^ = (V>i,j + V , i,*)/V2 ) 72, i = (ip2 
ipl i )/(v / 2i), we have 



(23) 
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— 2 J2 sin k 
2J2 sin k 



r 



Am 

V i>B,-k I 



The excited energies of the edge states are 



2J 2 sin k + J J{ a + J{ h - 2J x . a Ji.b cos k 



E- = 2J 2 sin k - J J( a + J( b - 2 J lia Ji,6 cos k 



If Ji.a — Ji.b (or A = B), the system has translational 
symmetry and the edge modes becomes gaplcss. In gen- 
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oral, due to J2 << Ji, a ,Ji,bi we have the full gapped 
edge modes. In Fig. 12, we set J 2 = 0.0125Ji. a and 
Ji. b = 0.5Ji, a . 

IV. CONCLUSION 

In this paper we study the symmetry protected Majo- 
rana edge states for the Z2 topological order of the Wcn- 
plaquette model and those of the toric-code model and 
calculate the dispersion of the Majorana edge states. We 
found that for the Wen-plaquette model the nodal points 
are fixed at k = and k = n on an edge which is pro- 
tected by the Z2 topological invariants and translational 
invariancc. Due to the interference effect between the 
Majorana fcrmions on two boundaries, the nodal points 



of the Wen-plaquctte model are not fixed at k = 0/n. 
For strong interference case, the edge states become full 
gapped. While, for the toric-code model with A 7^ B, 
there is no translational symmetry along the edges. Then 
the edge modes becomes gapped. If we recover the trans- 
lational symmetry by setting A = B, the edge states of 
the toric-code model have same properties to those of the 
Wen-plaqucttc model. 
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